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Abstract 

CDF has observed a top forward-backward asymmetry discrepant with the Standard Model pre- 
diction at 3.4cj. We analyze models that could generate the asymmetry, including flavor-violating 
T^'s, horizontal Z'jjS, triplet and sextet diquarks, and axigluons. We consider the detailed pre- 
dictions of these models for the invariant mass and rapidity distributions of the asymmetry at 
the parton level, comparing against the unfolded parton-Ievel CDF results. While all models can 
reproduce the asymmetry with the appropriate choice of mass and couplings, it appears at first ex- 
amination that the extracted parton-level invariant mass distribution for all models are in conflict 
with Tevatron observations. We show on closer examination, however, that tt events in Z'jj and 
W models have considerably lower selection efficiencies in high invariant mass bins as compared to 
the Standard Model, so that W' , Z'^, and axigluon models can generate the observed asymmetry 
while being consistent with the total cross-section and invariant mass spectrum. Triplet and sextet 
models have greater difficulty producing the observed asymmetry while remaining consistent with 
the total cross-section and invariant mass distribution. To more directly match the models and 
the CDF results, we proceed to decay and reconstruct the tops, comparing our results against the 
"raw" CDF asymmetry and invariant mass distributions. We find that the models that successfully 
generate the corrected CDF asymmetry at the parton level reproduce very well the more finely 
binned uncorrected asymmetry. Finally, we discuss the early LHC reach for discovery of these 
models, based on our previous analysis |48j . 
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I. INTRODUCTION 

Recently, there has been an apparent anomaly in the top sector: the observation by 
the CDF experiment of a top forward-backward asymmetry (Afb) [I], where the forward- 
backward asymmetry in a particular invariant mass bin, Mfj j, is defined by 

,u(^ . _ N{Ay > 0, Mu,) - N{Ay < 0, Mg,) , 
^ "'^^ N{Ay > 0, Mu,) + N{Ay < 0, M«-,) ' ^ ^ 

with Ay the rapidity diff'erence between a top and an anti-top. The recent CDF anlaysis 

shows Afb = 0.475 ± 0.114 for Mu > 450 GeV [Ij, while the Next-to-Leading Order (NLO) 

Standard Model (SM) predicts much lower values 0.088 ± 0.013 pHH], corresponding to a 

3.4cr deviation. The DO collaboration also observes a larger than predicted asymmetry [6]. 

Since the SM prediction for the top-pair production cross section is in relatively good 

agreement with observation, a new physics model must generate the large forward-backward 
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asymmetry without disturbing the total cross-section or observed invariant mass spectrum of 
tt production. In order to do this, many models assume an additional tree-level contribution 
from the exchange of a new particle in a way that maximizes the effect on the forward- 
backward asymmetry while minimizing the effect on the overall production cross section. 
Models proposed thus far in the literature that generate the observed Tevatron Aps at 
tree level while not grossly disrupting the top pair production cross section fall into two 
categories according to the nature of the new particle exchange: (i) s-channel exchange of 
vector mediators with axial couplings {e.g. axigluon models) or (ii) t-channel exchange 

of flavor- violating mediators pUHS^ . Comparative studies of these models have also been 
carried out |35] - H5] . The s-channel mediators tend to have maximally axial couplings, while 
the t-channel mediators tend to be maximally flavor violating, connecting a light quark to 
the top quark. Recently, it has been pointed out that such maximal flavor violation can also 
explain anomalies in the Bg and systems when the 6-quark is coupled as well ^61 H7] . 

The recent CDF analysis [T] greatly extends the experimental information about the 
asymmetry. Due to much improved statistics of top-pair production at the Tevatron, the 
analysis shows event distributions in more detail. In particular, CDF has now collected 
enough data to give the forward-backward asymmetry with respect to the invariant mass, 
Mu, and the rapidity, y, of reconstructed tops, and to compare these distributions for various 
data subsets, such as the 4- and 5-jet samples. This alows us to reassess the viability of 
previously suggested theoretical models that address the Aps anomaly. For example, in 
the high invariant mass bin, Mu > 450 GeV, the asymmetry is very large, ~ 50%, and 
hence some models that marginally explain the previous Aps will be challenged. Diquark 
models, for example, have great difficulty in producing a large enough asymmetry without 
generating an overly large correction to the total production cross-section. 

In addition, although the top forward-backward asymmetry and related variables are 
straight-forwardly defined observables, the interpretation of experimental observation to 
confirm or falsify models must be done with care due to the indirect nature of top quark 
identification. In spite of the effort of the CDF group to show an "unfolded" ApB for direct 
comparison with theory predictions, a broad coverage of models is needed to justify its 
model independence, especially where selection effects can come into play. Therefore, it is 
necessary to compare direct experimental signatures to expectations within a wide variety 
of models. 
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The purpose of this paper is to give a comprehensive analysis of a variety of models that 
have been proposed in the literature to explain the Aps anomaly in order to (1) reassess 
the viability of such models given the new experimental data and (2) investigate subtleties 
associated with unfolding the "raw" ApB to a parton level Aps for comparison to models.^ 
In order to do the latter, we must decay the tops, sending the events through a detector 
simulation, and reconstruct the tops utilizing a x^-based top reconstruction algorithm. Be- 
cause, unlike previous studies, we do fully reconstruct the tops and simulate detector effects, 
we are able to make comparisons to the raw CDF asymmetries. 

Comparing our results against the raw asymmetries and invariant mass spectra also 
enables us to assess the model dependence of event selection efficiencies. We find that 
these efficiencies can dramatically affect the extracted parton level total and differential 
cross-sections.^ For example, Z'^j and W models that produce large asymmetries have event 
selection efficiencies in the high invariant mass bins that are about half or less than what 
a Standard Model tt Monte Carlo would predict. This must be taken into account when 
comparing a model against parton level results. For example, while the parton level invariant 
mass spectra of Z'^j and W models appear badly in conflict with the reconstructed invariant 
mass spectra from CDF at high invariant mass, we flnd that once the selection efficiencies 
are taken into account, the agreement between the data and model is good. 

This paper is complementary to our earlier paper in which we examined the reach of 
the LHC at 7 TeV to discover the mediators of t-channel physics generating the Tevatron 
asymmetry. We found there that a top-flavor violating state light enough and with large 
enough couplings to generate the asymmetry will be rapidly discoverable at the LHC via a 
search for tj resonances in ttj events |18]. 

The outline of this paper is as follows. In the next section we summarize and discuss 
the qualitative features of the classes of models that have been shown to be capable of 
generating the Tevatron top ApB- In the following section, we flrst show our parton- level 
comparison of various models to the unfolded Tevatron results. Then we show detector-level 
comparison utilizing a reconstructed top sample. Next we present a lepton asymmetry from 

^ We use "raw" to refer to the measurement uncorrected for detector effects, but with background sub- 
tracted, as presented in the CDF paper. 

^ We use "parton level" as in the CDF Aps paper, where the term refers to de-convolving event selection 
efficiencies, detector efficiencies, jet algorithms, background, etc., from the underlying physics jlj. 
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(a) s-channel qq (b) s-channel gg (c) t-channel qq 

Figure 1: Tree level ti production diagram with mediator M exchange. 

fully leptonic ti events, which has recently been discussed in a CDF note Lastly, we 
discuss the LHC reach for discovering such states, based on the analysis of |48j . 

II. MODELS 

The Leading Order (LO) SM tree-level amplitude for tt production does not generate a 
forward-backward asymmetry. In the SM, a small positive top forward-backward asymmetry 
is generated through interference between a one-loop box diagram and a LO tree level 
diagram, Afb{Mu < 450GeV) = 0.040 ± 0.006, Afb{Mu > 450 GeV) = 0.088 ± O.OIS.^ 
Since the SM contribution is generated at NLO, if there is an additional LO tree-level 
contribution from new physics, it can easily dominate. 

Such LO diagrams are of the form of those in Fig. ([T]). They can be either s-channel 
(Fig. ([T^) and ^i)) or t-channel (Fig. s-channel mediators couple directly to light 

flavors and gluons, and therefore the mediator masses must be large enough to evade dijet 
resonance search constraints pTl [T7] . To maximize the contribution to Aps, such a model 
must have a big axial coupling. 

On the other hand, t-channel models should have large flavor violation between the light 
and the top generations, as can be seen in Fig. Large flavor violation is experimentally 
allowed even for low mass mediators, M, as long as new couplings between light generations 
and left-handed quarks is suppressed; then strong limits on flavor violation and from dijet 

^ Interference between initial state gluon radiation and final state gluon radiation makes a very small 
negative contribution to the asymmetry. 
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resonance searches are avoided. Additionally, the same-sign top signature search limit prefers 
M to be a non-self-conjugate state ^U\. Therefore, ordinary Z' models run into difficulty. 
Here, to avoid same-sign top constraints, we consider horizontal Z^s with flavor charge. 
Color exotic states and W's can also satisfy the requirement. 

In the following sections, we summarize the defining Lagrangian of t-channel W, Z'^j, 
triplet scalar, sextet scalar and s-channel axigluon models and present the tree-level differ- 
ential cross sections, ^^Mzt^_ 



A. Flavor-Changing W', Z' 



The Lagrangian for a flavor- violating Z' interaction is 



1 - 



C = —tYigiPL + 9rPr)uZI + h.c. 



giving rise to a scattering cross-section 

da /3 



{•Asm + Aint + At 



where 



dcos9 327rs 

Asm = ^{1 + 4 + -:^, 
9 s 



qJ ' 



with ce = (3 cos 6 and /3 = a/1 — Am^/s. The new physics contributions are 



A: 



int 



9 siz' 



2ul + 2sml + ^{il + sml) 



(2) 



(3) 



(4) 



(5) 



A. 



sq 



2i%, 



(91 + 9rM + 29i9RKs - 2mi) + ^{91 + 9itr{4' + 4smi, 



(6) 



with ti = t — mf and Ui = u — ml. The Mandelstam variables are related to the scattering 
angles via i = — s(l — ce)/2 + and u = — s(l + C6))/2 + m^. Note that the Lagrangian has 
been deflned with a \/2 with respect to some other conventions in the literature. Similar 
expressions hold for the flavor-violating W via the interaction Lagrangian 



C = -^d^igLPL + 9RPR)tW'^ + h.c. 



(7) 
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B. Color Triplets and Sextets 



The quantum numbers of the color triplet and sextet are 



(3, 1)4/3 (6, 1)4/3, 



and their interactions with up and top quarks is given by 



£^ = rt'VigiPL + gRPR)u. 



This gives rise to a scattering cross-section 
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9 ul 



(8) 



(9) 



(10) 



where C(o) = 1(— 1) for triplets (sextets) [22l is a color factor that comes from the 
interference of new t-channel physics with the s-channel gluon. The color factor C(2) comes 
from the squared new t-channel physics term and is equal to C(2) = 3/2 for sextets and 
C(2) = 3/4 for triplets. We have also defined g = a/ {g1 + q'r)/'^- 

C. Color Octet 



The exotic gluon couples to light quarks through 

Ca.^ = gs {qT^YiglPL + gWq + tT^r{g\PL + g'^Pn)t) G'^. 



(11) 



Note the inclusion of the QCD coupling constant, gs, in the interaction. The scattering 
cross-sections calculated through these interactions are 



A.. 



9: 



sis — m 



G' 



int 



9 {s- ml,Y + ml,Tl, 



igl + glWL + gn) 



(o /Q2^ I p idL 9r){9l 9r) ^ 2 
^ ^^^^ i9l + 9WL + 9^n) '^'r 



(12) 



A. 



.1^ 



+ 9'n'){9i' + 9'n 



sq 



m 



t 2^ 



G'-^ G' 
-.t „t 



2^ ^9l9R n (9f-9'R){9i 

2 , ^ 2 + ^ / a2 



.'?2\/ t 2 



l + (l-/3 



t 2\ 

9r ] 



9'l+9'r ■ '{9r + 9'RWL+9'R 



C0 + Cg 



(13) 



As does CDF, we consider the case where the couplings of the vector color octets are 
purely axial, so gy = {g\ + 5'^)/2 = and gy = {g\ + g^ji)/2 = 0, and where the axial 



coupling of the boson to light quarks is positive and opposite the coupling of the boson to 
tops, so g\ = {gj^ — gi)/2 = —g\ = {—gfi + 5'l)/2. This axigluon case leads to the largest 
positive contribution to the asymmetry per contribution to the total cross-section. 

III. PARTON LEVEL TEVATRON TOP FORWARD-BACKWARD ASYMME- 
TRY 

We now simulate the models described in the previous section, and we will use the 
formulae presented there to discuss the qualitative features. 

We begin by analyzing models for top asymmetry generation at the parton level. "Parton 
level" as used in the CDF ApB paper refers to de-convolving event selection efficiencies, 
detector efficiencies, jet algorithms, background, etc., from the underlying physics [IJ. Thus 
for the parton level analysis, we simulate and compare our results after showering in PYTHIA 
but before folding in detector effects. To do this, we use MadGraph/MadEvent with matrix 
element / parton shower (ME/PS) matching in the MLM scheme, which was implemented 
in the MadGraph/MadEvent package using PYTHIA. The events were generated using a fixed 
renormalization scale and factorization scale of 200 GeV.^ MadGraphS vO . 6 . 1 / MadEvent 
4 . 4 . 44 with QCUT = 30 and xqcut = 20 was used to generate all signal and standard model 
events. Showering and matching are done in order to improve accuracy, including the effects 
of single mediator production in inclusive tt events. 

A. Cross-sections and Efficiencies 

We make a few observations about total cross-sections and invariant mass distributions 
before delving into a detailed analysis of the asymmetry. Fig. ^ shows the parton level 
A^p^ for events with rritt > 450 GeV versus total leading order pp ^ ti + or 1 jet 
cross-section for a swath of flavor-changing W, Z'j^, triplet, sextet models, and axigluon 
models. A comprehensive scan of models was carried out, and only representative points 

^ The current version of MadGraph uses as{mz) = 0.13, which is substantiahy larger than the current 
measured value of as = 0.118. Therefore our choice of the renormalization scale is effectively lower than 
the nominal scale, assuming as(™z) = 0.118. Our choice reproduces well the known theoretical LO 
Standard Model value for the tt cross-section at Tevatron. 
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Figure 2: Parton level forward-backward asymmetry for events with invariant mass greater than 
450 GeV versus leading-order cross-section in picobarns for various models. The mass (in GeV) and 
coupling of the mediator are indicated by (co^pUng) • ^ comprehensive scan of models was carried 
out, and only representative points are shown. The coupling shown for axigluons is g*^ = —g\ 
(supposing gy = gy = 0) while the coupling shown for Z'^ and W models is gn (assuming = 0) 



and the quoted couplings for triplets and sextets is g = \J {g'l + For comparison against 

observations, the horizontal shaded area lies in the itlcr region of the measured (parton-level) value 
of ApsiMfi > 450 GeV). The vertical lines lie at the central value of the CDF tt production cross- 
section (7.5 lb 0.48 pb), divided by a K-factor of 1.3 or 1. The SM marker lies at the value of the 
LO Standard Model cross-section and the NLO value for the SM forward-backward asymmetry. 
Note that care must be taken when comparing the new physics cross-sections against the Standard 
Model cross-section, as the selection efficiencies for the new physics models can be lower. This is 
discussed in more detail in the main text. 

are shown. The horizontal shaded band lies at ±la values of the observed asymmetry 
in the high invariant mass bin, > 450 GeV. The vertical dashed lines correspond to 
the combined ti cross-section from CDF with 4.6 fb^^; CDF measures a cross-section of 
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2TcVAiigluoii,g/=-j„'=2.4 
,00GeVZ„',g=3.4 



Figure 3: Efficiencies after showering events with PYTHIA and making the cuts detailed in section IV 
on final state leptons and jets, with the hatched regions corresponding to la errors based on the 
limited statistics of the sample. No detector effects have been taken into account here. Black bars 
are for LO Standard Model. All samples are LO matched. 

7.5 ± 0.31(stat) ± 0.34(syst) ± 0.15(lumi), assuming a top mass of nit = 172.5 GeV [50] . 
The predicted next-to-leading-order (NLO) SM cross-section at the value of the top mass 
we assumed in simulations, rrit = 174.3, is about 7.2 pb [51j, whereas we find the LO SM 
cross-section is 5.6 pb, implying a SM K-factor of about 7.2/5.6 ~ 1.3. Of course, the NLO 
corrections to the new physics have not been calculated, so any comparison between the 
observed cross-section and the ti production cross-section is subject to some uncertainty. We 
do choose to show in Fig. ([2]), however, the central value of the combined CDF ti production 
cross-section (7.5 pb) divided by the SM K-factor when comparing to the leading-order (LO) 
tt production cross-section of SM plus new physics against the observed production cross- 
section. From the figure it is clear that, in general, excepting the Z'^j and axigluon models, 
models with couplings that are small enough to be in accord with the observed cross-section 
do not produce a large enough asymmetry in the high ti invariant mass bin. 

This statement requires a strong qualification, however, which we investigate in detail be- 
low after carrying out a reconstruction of the top samples in these models. The qualification 
is that the efficiency for a ti event to pass cuts (the same as those used in the CDF analysis 



and our detector level analysis below, see Sec. (IV)) is strongly model dependent for cases 



where there is a large asymmetry. This is shown in Fig. ([3]), where we see that the efficiencies 
to pass cuts (after showering and jet clustering but no detector simulation) is suppressed by 
more than a factor of two relative to the Standard Model for the 400 GeV model shown 
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LQ Standarti Model 




and by about a factor of 1.5 for the 800 GeV Z'^ and the W model. The reason for this 
becomes clear after examining the distribution of events in top pseudorapidity, r/^, and and 
tt invariant mass, m^f, as shown in Fig. Q. In order to generate a large asymmetry, the rft 
distribution must be skewed significantly with respect to the SM distribution. To generate a 
very large asymmetry in the high invariant mass bin, the distribution must be more skewed 
at high invariant mass. The distributions for the Z'^ and W models are so skewed at high 
invariant mass that the peak of the distribution lies close to the rjt = 1 line. Thus a cut on 
lepton pseudorapidity of \ri\ < 1 and jet rapidity of |?7| < 2 ends up cutting out a signifi- 
cantly greater fraction of events at high invariant mass than in the SM case. Importantly, in 
unfolding the differential ti cross-section, assumptions about event selection efficiencies must 
be made; the assumption is that actual event selection efficiencies do not differ substantially 
from the Standard Model efficiencies [52]. [53] . This assumption clearly breaks down in the 
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rritt < 450 GeV 


nitt > 450 GeV 




Ay < 


Ay > 


Ay < 


Ay > 


Model 


eff. 

JJ 


r 


eff. 

JJ 


r 


eff. 

JJ 


r 


eff. 

JJ 


r 


SM 


0.079 


0.31 


0.078 


0.3 


0.092 


0.2 


0.089 


0.19 




0.078 


0.21 


0.076 


0.22 


0.088 


0.15 


0.063 


0.42 


W 


0.079 


0.23 


0.077 


0.27 


0.095 


0.16 


0.075 


0.34 


Triplet 


0.084 


0.18 


0.083 


0.23 


0.103 


0.2 


0.095 


0.39 


Sextet 


0.075 


0.26 


0.073 


0.28 


0.087 


0.19 


0.08 


0.27 


Axigluon 


0.079 


0.26 


0.077 


0.31 


0.096 


0.14 


0.086 


0.28 


r7l 2 


0.074 


0.18 


0.072 


0.19 


0.089 


0.16 


0.069 


0.47 



Table I: Parton level efficiencies, eff., and bin fractions, r, for the Standard Model (SM) and 
for benchmark models Z'j^^ ( 400 GeV, gn = 1.75), Z'^^^ ( 800 GeV, gn = 3.4), W (400 GeV, 
gn = 2.55), Triplet (600 GeV, g = 4.4), Sextet (1.4 TeV, g = 4.0), and Axigluon (2 TeV, g\ = 

_ t =2 4) Here eff = # events in bin after cuts i ^ _ # events in bin after cuts 
"A ' )' JJ' events in bin before cuts total # events after cuts 

case of the Z'^ and W . We investigate this effect in the reconstructed sample in the next 
section. For now, we compare the parton level asymmetries and cross-sections, and note 
caveats where this effect is important. 

Efficiencies can also affect the forward-backward asymmetry. If, for example, the ef- 
ficiency for an event to pass cuts is lower for events with Ay > than for events with 
Ay < 0, then cuts will wash out the asymmetry. We show the efficiency to pass cuts (after 
showering and jet clustering but no detector simulation) in four mjf. Ay bins for the SM 
and for several benchmark models in Table [H The difference between SM efficiencies and 
new physics model efficiencies is not as great for the coarse uiti binning, implying that the 
effect on the coarse binned asymmetries reported by CDF will not be as great as for the 
more finely binned invariant mass distribution. However, we also see from the table that 
the efficiencies for Ay > are smaller for some new models, such as the Z'^ and W than 
for the SM and the other new models, implying more washout of the asymmetry for the 
and W , which should be compensated. These effects will be seen when we compare our 
parton level asymmetries to the reconstructed asymmetries. The overall point here is that 
the efficiencies have some effect on unfolding the parton level asymmetries, but they are not 
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Model 


Mass(GeV), coupling, cross-section(pb) 


FV W 


200, 1.4, 7.1 


300, 1.8, 6.5 


400, 2.4, 6.9 


400, 2.6, 7.7 


600, 3.4, 6.9 


600, 3.6, 7.4 


FV Z'jj 


300, 1.4, 6. 


400, 1.6, 5.1 


400, 1.8, 6.2 


600, 2.4, 5.6 


800, 3.2, 6. 


800, 3.4, 6.8 


Triplet 


400, 3., 7.9 


400, 3.2, 9.5 


600, 3.6, 6.7 


600, 3.8, 7.4 


600, 4., 8.4 




Sextet 


600, 2., 8.1 


800, 2.4, 7.8 


1000, 3., 8. 


1200, 3., 7.1 


1200, 3.4, 7.7 


1400, 4., 7.7 


Axigluon 


2000, 2., 5.7 


2000, 2.4, 5.8 


2000, 3.2, 6. 


2200, 3.2, 5.8 







Table II: Summary of benchmark models. The coupling shown for axigluons is (7^ = —g\ (supposing 
9v = Qv = 0) while the coupling shown for Z'^ and W' models is qr (assuming ql = 0) and the 
quoted couplings for triplets and sextets \s g = Y^(ff|^^r5^)/2. For each model we consider, we 
include mass, coupling and total leading order matched + or 1 jet production cross-section as 
calculated with MadGraph/MadEvent/Pythia. The cross-sections above can be compared to the 
cross-section we obtained for the Standard Model using the same cuts and SM input parameters: 
5.6 pb. No K-factors have been included in these quoted cross-sections. Note that care must be 
taken when comparing the new physics cross-sections against the Standard Model cross-section, as 
the selection efficiencies for the new physics models can be lower. This is discussed in more detail 
in the main text. 

as large an effect as on the invariant mass distribution. 
B. Parton Level Asymmetries 

We now show the parton level asymmetries for each of the benchmark models. We choose 
several benchmark masses/couplings for Z^, W , triplet, sextet, and axigluon models that 
give rise to large forward-backward asymmetries without generating too large a total cross- 
section. These benchmark models are listed in Table [Tll 

We show in Fig. ([s]) the forward-backward asymmetry in Mti and Ay, comparing the 
reconstructed parton level asymmetry of the Tevatron to our simulated matched sample. 
The hatched regions correspond to the la errors based on the limited statistics of the 
sample. We can see that the horizontal Z' and W can give good fits in both the high and 
low invariant mass bins, and in low and high Ay. At least for the two bin case, the models 
also appear to give the correct shape as a function of Mti and Ay. 
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^^^^^^ 



200 GeV,g=l.4 
300 GeV,g=1.8 
400 GeV, g=2.4 
400 GeV, g=2.6 
600 GeV, g=3.4 
600 GeV, g=3.6 



200GeV, g=1.4 
300GeV, g=1.8 
400 GeV, g=2,4 
400 GeV, g=2,6 
600 GeV, g=3,4 
600 GeV, g=3.6 



Figure 5: ApB for Z^, W' models with couplings as indicated in the legend (with g = gji, gi = 0), 
with the hatched regions corresponding to la errors based on the limited statistics of the sample. 
The contribution to Aps includes both t-channel Z^, W exchange and single Z'^, W' production. 
The red bars are the CDF observation with la errors, while the blue bars indicate the NLO SM 
contribution from [1] , which has not been included in the LO contribution calculated via MadGraph 
and PYTHIA. The last bin includes all events with m^j > 450 GeV and |Ay| > 1, respectively. 

A similar analysis is carried out for triplets, sextets and axigluons in Fig. (|6|. While 
triplet and sextet models can marginally reproduce the asymmetries at LO, the rise between 
the low and high invariant mass bins and low and high rapidity bins is not as pronounced for 
the triplets and sextets as for the W and Z'^. The reason for this in the sextet and triplet 
cases can, for example, be easily extracted from the analytical expressions, Eqs. ([s]), (|6]), ( 10 ). 



At high invariant mass, the scattering amplitude is dominated by the squared term. There 



Ut — in Eq. (10), and the effect of the Agq on the asymmetry vanishes. By contrast, the 
Z'^, W are dominated by u^/ij^f ~ (1 + cg)^/ (1 — cg)"^ which retains a contribution to the 
asymmetry at high invariant mass. 

The axigluon models tend to significantly underproduce the asymmetry in the high in- 
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Triplet 



Triplet 



iHrHrHrHvH:;] GeV, g=3. 

^^^^^ 400GeV,g=3.2 

'(c^i/^^^SZSI^^^^Si fiOO GeV, g=3.6 

^^^^^^^^^^^^^^^^ 600GeV,g=3.8 

: ;V: : ;^ 600 GeV, g=4. 
^7X7 7;7v; 7^77; 7 >v>^-^^^^^7^ 



Sextet 



600 GeV, g=2. 
^^^^^^ 800GeV,g=2.4 
^^J^iX^^l^i! 1000 GeV, g=3. 
Y//////////////A 1200 GeV, g=3. 

: ; : ; :V: ;V^^ 1200 GeV, g=3.4 
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Figure 6: ^4^?^ for triplet, sextet and axigluon models with couplings as indicated in the legend, 
with the hatched regions corresponding to \a errors based on the limited statistics of the sample. 
The red bars are the CDF observation with \a errors, while the blue bars indicate the NLO SM 
contribution, which has not been included in the LO contribution calculated via MadGraph and 
PYTHIA. The last bin includes all events with m^f > 450 GeV and |Ay| > 1, respectively. 

variant mass window. Choosing a larger coupling does not give rise to a larger asymmetry in 
the high invariant mass bin because of width effects, and the axigluon mass cannot be low- 
ered in order to compensate on account of dijet constraints [T7]. Thus we see that axigluon 
models have greater difficulty than and W for reproducing the observations. Some of 
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Figure 7: -^^^ for the benchmark models appearing in Figs. (5, 6). The Tevatron measured cross 
section (red crosses, from [52J), and LO SM cross-section with the same SM parameters and fixed 
renormahzation scale used to generate benchmark model events are also shown. No K-factors are 
applied. 



these constraints can be relaxed somewhat by moving away from the point g\ = —g\ [TT] - 
In addition, on account of the large couplings present in these models, NLO corrections to 
the new physics must be considered in order to draw firm conclusions. 

Before moving on to the fully reconstructed sample, we compare the invariant mass 
distributions of the LO PYTHIA results against the observations in Fig. ([T]). No K-factors 
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for NLO corrections have been applied. All the models overproduce the extracted invariant 
mass spectrum in the high mass bins. Here again, however, the caveat must be applied that 
the Z'^ and W models have lower selection efficiencies in the high invariant mass bins, so 
that one expects the discrepancy in the bins above rriti ^ 500 GeV to be greatly reduced 
for these models. On the other hand, the sextets and triplets severely overproduce the 
observed number of events and, based on Fig. Q, are not helped by having lower selection 
efficiency in the high invariant mass bins. The axigluon models do not have as severe an 
overproduction problem, but also do not generate a large asymmetry, as can be seen in 
Fig. ([6]). As commented earlier, it is of course possible that NLO corrections from the 
new physics will lead to significant changes in these distributions on account of the large 
couplings. 

We now turn to the reconstructed sample from which it will be possible to make more 
quantitative statements about the observed versus model-dependent predicted invariant 
mass spectra and asymmetries. 

IV. FULLY RECONSTRUCTED ASYMMETRY AND INVARIANT MASS DIS- 
TRIBUTIONS 

To reconstruct the invariant mass spectrum and check the model dependence of the tt 
parton-level asymmetry extracted in [1] for the class of models discussed here, we can send 
the showered tt events through PGS. To select the tt signal, we take the same requirements 
as CDF in their analysis: 

• Exactly one electron or muon with > 20 GeV and |?7| < 1.0. 

• Photon and r veto. 

• At least four jets with > 20 GeV and \r\\ < 2.0, with at least one of the jets having 
a 6-tag. 

• E^'"' > 20 GeV. 

We must then fully reconstruct the decayed tops. We do a likelihood analysis on the lepton 
and jet kinematics to the ti hypothesis, using the algorithm described in our previous paper 
[is] . The top is reconstructed out of the four hardest jets in the event. In order to gain 
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enough statistics in the high invariant mass bin to rehably compare our results against the 
reconstructed CDF asymmetry, we generate 5 miUion ti events per model. Approximately 
2% of these events survive the cuts. In contrast to our previous analysis [48J, but in ac- 
cordance with the CDF analysis, we place no cut on the reconstruction of the tops. We 
explore later the effect of the cut on the size of the asymmetry. 

Because of the large numbers of simulated and reconstructed events required, we consider 
only a representative subset of the models analyzed at the parton level in the previous 
section. We choose a 400 GeV Z'^ with qr = 1.75, a 400 GeV W with qr = 2.55, a 600 
GeV triplet with g = 4.0, a 1.4 TeV sextet with g = 4.0 and a 2.0 TeV axigluon with 
g\ = —g\ = 2.4.^ As we will see, the CDF extraction of the parton level asymmetries and 
invariant mass spectra is somewhat model dependent, on account of the model dependent 
efficiencies shown in Fig. ^ and Table (|l]), as well as detector effects. 

We begin by comparing our reconstructed results against the CDF results for several 
models in Fig. The shown there is as defined in Eq. ([I]), with the top and anti-top 
identified by the sign of the lepton. We see that models that reproduce the parton level 
asymmetry match well against the fully reconstructed CDF asymmetry. The Z'^j and W 
models, however, receive a larger upward correction upon unfolding to the parton level than 
would be expected using SM efficiencies. The reason for this is clear from Table (|l]): the 
efficiencies in the high invariant mass bin with Ay > for and W models are lower 
than for the SM, leading to a greater washout of the asymmetry at the detector level. We 
also observe that the axigluon, while appearing to reproduce the reconstructed asymmetry 
marginally, tends to underproduce the unfolded parton level asymmetry as seen in Fig. (|6]). 
We can also look at the partitioned asymmetry defined by 

,tt( ^ ^ N{{y, - yn) > 0, q, Mg,) - N{{y, - y^) < 0, g, Mg,) 

""^ Ar((y,-y,)>0,g,M,,-,) + Ar((y,-y,)<0,g,M„-,)' ^ > 

where ye, is the rapidity of the leptonic top, yu is the rapidity of the hadronic top, and q is 
the charge of the lepton. The asymmetries obtained in this way are shown in Fig. ([9]). 



^ Note that the W' and Z'^ models wiU require a triplet or higher Higgs representation in order to evade 
dijet constraints for the flavor-conserving Z's, which exist in these models. On the other hand, we also 
compare a heavier 800 GeV Z'^ to the lighter W' and Z^, and find that it overproduces the high invariant 
mass spectrum. 
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Figure 8: A^p^iMiii) for W' , Z'^, triplet, sextet and axigluon models. Red crosses are the CDF 
values reconstructed from data. Blue crosses are the MC@NLO expectation. The last bin includes 
all events with m^f > 700 GeV. 
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for W\ Z'^, triplet, sextet, and axigluon models. Here, 
the data is divided according to the charge, of the lepton in the event. Black {q > 0) and red 
{q < 0) crosses are the CDF values reconstructed from data (with background subtracted). 



Figure 9: Ap^{q, M^) as defined in Eq.( 14 
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Figure 10: Number of expected tt events with 5.3 fb~^ at the Tevatron, distributed over m^j. 
Events were passed through PGS and then tops were reconstructed using the algorithm detailed 
in [lEj. The red bars indicate CDF's measurement, with expected background (as estimated by 
CDF) subtracted. The green histogram is our SM sample and the purple histograms represent 
model samples. For SM and model samples, a fixed number of events were generated, and then 
event counts were scaled appropriately for 5.3 fb^^ integrated luminosity. 

We also compare the reconstructed invariant mass spectrum against that reported in [1]. 
Since the efficiencies are lower for the Z'j^ and W than for the Standard Model, we may 
expect these models to agree better with the observations than suggested by the parton 
level invariant mass spectra shown in Fig. ([T]). We compare in Fig. (10) the simulated recon- 



structed invariant mass spectrum against that reported in p?| . First we note the discrepancy 
at low invariant mass between all models (including the SM) and the observations, which 
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Figure 11: ^4^^ and ^4.^^ in low (m^j < 450 GeV) and high (m^j > 450 GeV) ti invariant mass 
bins for Z^, VF', triplet, sextet, and axigluon models. Red / Blue crosses are the CDF values 
reconstructed from data in the lab / CM frames. Purple / Green crosses indicate the SM NLO 
predictions. 

we attribute to NLO corrections and to a difference between the PGS detector simulation 
and the CDF simulation. However, we can see the effects of the efficiencies in the high 
invariant mass bins noted in Fig. ([s]). For example, we can see the efficiency correction does 
seem to bring the W model into agreement with the SM. On the other hand, the triplet 
model largely and almost uniformly overproduces the invariant mass distribution in all bins. 
While the sextet model appears in better agreement at high invariant mass, it underpro- 
duces the observed asymmetry as shown in Fig. ([s]). In general, triplet and sextet models 
have greater difficulty producing the observed asymmetry while remaining consistent with 
the total cross-section and invariant mass distribution, as emphasized by Figs. ([2||3]). 

There are two other comparisons that we are able to do with our fully reconstructed 
asymmetry. We are able to compare the center-of-mass versus lab frame asymmetries, which 



is shown in Fig. (11). While the models give rise to some difference between the CM and 



lab frames, the difference is less pronounced than what CDF observes. We can also compare 



the asymmetries in the four and five jet samples, as shown in Fig. (12). Here we see some 
washout of the asymmetry in the 5 jet sample, an effect that is observed in the CDF data. 

Lastly, though this is not considered in detail in the CDF analysis (a value for the raw 
asymmetry after a cut of 3 is presented in Table XIV of [Ij), it is interesting to observe 
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Figure 12: vl^^ in low (m^j < 450 GeV) and high (m^j > 450 GeV) invariant mass bins for 
benchmark models. The simulated data samples were partitioned according to whether the event 
had more than five jets with > 20 GeV and |r/| < 2. The CDF measured lab frame asymmetry 
for 4-jet and 5+-jet samples is shown as red crosses and black crosses, respectively. 
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Figure 13: J^p^{Mii^) for Z'^ model with various cuts on the it reconstruction x^. The last bin 
includes all events with mii > 700 GeV. Red crosses are the CDF values reconstructed from data. 
CDF used a likelihood algorithm for top reconstruction, but made no cut. 

the effect of a cut in the top reconstruction on the size of the asymmetry.^ We can see in 



Fig. ( 13 ) that a moderate cut increases the asymmetry especially in high invariant mass 
bins. 



^ Recall that we reconstruct tops by doing a fit on the lepton and jet kinematics to the tt hypothesis. 
The fit has three degrees of freedom. 
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Figure 14: Parton level lepton forward-backward asymmetry for tt events in which both tops decay 
leptonically, and in which both leptons pass the rapidity cuts corresponding to those used in the 
recent CDF analysis [49j. Red points are extracted from the results in [l9]. (a): ^_fb(Aj/j,+^_) ^ 
a function of |Ay| = - y^-]. (b): ApB{i^y^^^_^ 
i'^t' invariant mass. 
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V. LEPTON ASYMMETRY 



In order to avoid potential issues with the top reconstruction, one can also look at the 
asymmetries in di-leptons, where both tops decay leptonically. The CDF collaboration re- 
cently reported results from an analysis of di-leptonic ti events. In addition to reporting 
asymmetries obtained after reconstructing tops in events, they report the raw lepton asym- 
metry [l9]. We compare the raw lepton asymmetry in benchmark models at the showered 



parton level to the CDF measured value in Fig. (145,). To better compare with the CDF 



measurement, we show the asymmetry for only events that have electron rapidities in the 
range \7]\ < 1.1 or 1.2 < \ri\ < 2.8 and muons with rapidities in the range |?7| < 1, corre- 
sponding to the rapidity cuts placed on the leptons in their analysis. We also point out 
that examining the lepton asymmetry as a function of lepton-lepton invariant mass could be 
instructive. We have shown the lepton forward-backward asymmetry as a function of me+e- 
in Fig. (jlib). 
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VI. EARLY LHC REACH 



We now consider the feasibility of discovering such models in the early run of the LHC at 
7 TeV. To this end, we make use of our previous results [48J. According to these results, a 
200 GeV W with coupling of 1 should be discoverable at 3a with 1 fb^^ of data. For a 400 
GeV W, the coupling must be larger than 1.3, and for a 600 GeV W, the requirement is a 
coupling of 1.8. Thus we see all of the W models giving rise to large asymmetries should 
be observable with an fb~^ of data. 

Likewise, for the Z'^j model, a coupling larger than 0.7 is required to discover a 200 GeV 
state at 3 cr with 1 fb~^, while a coupling of 0.8 is required for a 400 GeV state, and a 
coupling of 1.2 for a 600 GeV state. For the triplets, a coupling of ~ 0.9 is required for a 3 
a discovery with 1 fb~^ for 400 GeV or lower masses; for a 600 GeV triplet, the requirement 
strengthens to requiring a coupling of 1.3. Similar types of constraints can be obtained for 
the sextet models. 

The broad conclusion here is that all of the t-channel models that we considered here 
to fit the Tevatron top forward-backward asymmetry should give rise to 3a excesses with 1 
fb^^ at the LHC in the context of a top-jet resonance search. 

According to the analysis in [17J, the axigluon benchmark models presented in this paper 
will be rapidly discoverable at the LHC through dijet events. 

VII. CONCLUSIONS 

We examined models of new physics that could generate the top forward-backward asym- 
metry. We considered W, Z'^ triplet and sextet di quarks, as well as axigluon models. We 
compared the asymmetries produced by these models to those observed at the Tevatron, and 
concluded that of the models that generate a large enough asymmetry in the invariant mass 
bin Mti > 450 GeV, the Z'^j, W and axigluon models are the only ones that do not hugely 
overproduce the total ti production cross-section. To bring the W models into agreement 
with the total ti production cross-section extracted at the Tevatron, we noted an impor- 
tant effect: the efficiency to select tt events from W models is significantly lower than for 
the Standard Model. This same effect is also helpful in improving the agreement between 
the invariant mass spectra of the W and models with the Standard Model predictions 
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(which agree with observations). Our result also differs from earlier studies on the diquark 
models which found that they could adequately produce the asymmetry without producing 
unduly large cross-sections. In order to further investigate the model-dependence in the ex- 
tracted asymmetry and invariant mass spectra, we then proceeded to decay the top quarks 
and simulate detector effects, reconstructing the tops via a likelihood based algorithm. This 
allowed us to compare our results against the raw CDF results in the asymmetry as well 
as invariant mass spectra. We found that when this was done, some W and Z'j^ models 
adequately reproduced the invariant mass spectra. It also allowed us to compare our results 
against the CDF results for lab versus center-of-mass frames, as well as the 4 jet versus 5 
jet asymmetries. We conclude that while the models reproduce the observed decrease in 
the asymmetry in the 5 jet sample, no appreciable difference occurs between the lab and 
center-of-mass frames. Lastly, we note that an LHC search at 7 TeV for top jet resonances 
could exclude at the > 3a level any Z'^j, W or diquark model that produces the Tevatron 
asymmetry. 

Acknowledgments: We thank Dan Amidei for many helpful conversations. 

Note added: While this work was being finalized, [Mj appeared, which explores sextet 
and triplet models. While our results agree quantitatively with theirs, our conclusions on 
the viability of these models for explaining the asymmetry are more pessimistic, on account 
of the total tt cross-section and invariant mass distribution. 
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